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Combining direct computations with invariance arguments, Taylor's constitutive equation for an 
emulsion can be extrapolated to high shear rates. We show that the resulting expression is con- 
sistent with the rigorous limits of small drop deformation and that it bears a strong similarity 
to an a priori unrelated rheological quantity, namely the dynamic (frequency dependent) linear 
shear response. More precisely, within a large parameter region the nonlinear steady-state shear 
viscosity is obtained from the real part of the complex dynamic viscosity, while the first normal 
stress difference is obtained from its imaginary part. Our experiments with a droplet phase of 
a binary polymer solution (alginate/caseinate) can be interpreted by an emulsion analogy. They 
indicate that the predicted similarity rule generalizes to the case of moderately viscoelastic con- 
stituents that obey the Cox-Merz rule. 



1. Introduction 

Apart from their technological importance, emulsions have served as model systems accessi- 
ble to rigorous theoretical modeling. The study of emulsions consisting of droplets of a liquid 
dispersed in another liquid has thus contributed substantially to our understanding of the rhe- 
ology of complex fluids. However, although major theoretical achievements date back to the 
beginning of the 20* century, further progress turned out to be difficult. The macroscopic rheo- 
logical properties of emulsions are determined by the reaction of the individual drops to the flow 
field, which in turn is modified by the presence of other drops. The mutual hydrodynamic in- 
teractions of drops complicates substantially the mathematical description. Moreover, depending 
on system parameters and flow type, droplets may break under steady flow conditions if a certain 
critical strain rate is exceeded. Rigorous calculations of the constitutive equation have therefore 
concentrated on very dilute emulsions and on conditions where drops are only weakly deformed. 
Sometimes, however, it is desirable to have an approximate expression, which — though not 
rigorous — can serve for practical purposes as a quantitative description in the parameter region 
beyond the ideal limits. As far as the dependence of the viscoelastic properties of an emulsion 
on the volume fraction (|) of the dispersed phase is concerned, such an approximation has been 



given by Oldroyd (1953 ). It is not rigorous beyond first order in iSf but serves well some practical 
purposes even at rather high volume fractions. It seems not reasonable to look for a comparably 
simple approximation for the dependence of shear viscosity r| on shear rate y that covers the 
whole range of parameters, where all kinds of difficult break-up scenarios are known to occur. 
In the next section, we propose instead a less predictive expression which contains an average 
drop size R (that may change with shear rate) as a phenomenological parameter The latter has 
to be determined independently either from theory or experiment. It turns out, however, that for 
a substantial range of viscosity ratios and shear rates, the expression for ri(Y) is to a large extent 
independent of morphology. For conditions, where drops do not break outside this region, we 
point out a similarity relation between this expression and the frequency dependent viscoelastic 
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moduli G'(a)), G"(co), similar to the Cox-Merz rule in polymer physics. More precisely, we show 
that in the limit of small drop deformation, the constitutive equation of an emulsion composed 
of Newtonian constituents of equal density can be obtained from the frequency dependent linear 
response to leading order in capillary number C and (reciprocal) viscosity ratio A-^ ' . And we 
argue that the identification is likely to represent a good approximation beyond this limit in a 
larger part of the C — A- parameter plane. Experimentally, this similarity relation can be tested 
directly, without interference of theoretical modeling, by comparing two independent sets of rhe- 
ological data. In summary, our theoretical discussion provokes two maj or em pirical questions. 
(1) If drops break: does the expression for the viscosity derived in Eq.( 2.18 ) describe the data 
with R the average drop size at a given shear rate? (2) If drops do not break below a certain char- 
acteristic capillary number C* (k): does the proposed similarity rule hold? To what extent does it 
generalize to non-Newtonian constituents? In Section |] we address mainly the second question 
by experiments with a quasi-static droplet phase of a mixture of moderately viscoelastic polymer 
solutions. 



2. Theory 

2. 1 . Taylor's constitutive equation for emulsions 

A common way to characterize the rheological properties of complex fluids such as emulsions, 
suspensions, and polymer solutions, is by means of a constitutive equation or an equation of 
state that relates the components pij + p5,j of the stress tensor to the rate-of-strain tensor eij. 
This relation can account for all the internal heterogeneity and the complexity and interactions of 
the constituents if only the system may be represented as a homogeneous fluid on macroscopic 
scales. The form of possible constitutive equations is restricted by general symmetry arguments. 



which provide guidelines for the construction of phenomenological expressions (Oldroyd 195C, 



1958). On the other hand, for special model systems the constitutive equations may be calculated 
directly at least for some restricted range of parameters. An early example for a direct computa- 
tion of the constitutive equation of a complex fluid is Einstein's formula 

il=il,(^l + ^c^^ (2.1) 

for the shear viscosity r| = pn/en of a dilute suspension (particle volume fraction (|) ^ 1). It is 
obtained by solving Stokes' equation for an infinite homogeneous fluid of viscosity ri^. containing 
a single solid sphere. For a sufficiently dilute suspension, the contributions of different particles 
to the overall viscosity r| can be added independently, giving an effect proportional to (|). In close 



analogy Taylor (1932) calculated r| for a steadily sheared dilute suspension of droplets of an 
incompressible liquid of viscosity r[ij = Xr[c in another incompressible liquid of viscosity rjc. For 
weakly deformed drops he obtained 

(2.2) 

which we abbreviate by rij- in the following. This expression includes Einstein's result as the 
limiting case of a highly viscous droplet, X ^ °°, As in Einstein's calculation, interactions of 
the drops are neglected. The result is independent of surface tension o, shear rate y, and drop 
radius R; i.e., it is a mere consequence of the presence of a certain amount (|) of dispersed drops, 
regardless of drop size and deformation (as long as the latter is small). Moreo ver, the dynami c 



(frequency dependent) linear re spon se of an emulsion has been calculated by pldroyd (1953[ ) 



His results are quoted in section 2.4 below. 



Under steady flow conditions, drop deformation itself is proportional to the magnitude of the 
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rate-of-strain tensor e,/. More precisely, for simple shear flow with constant shear rate y, the 
characteristic measure of drop deformation for given "k is the capillary number 



C 



(2.3) 



also introduced by Taylor (1934). It appears as dimensionless expansion parameter in aperturba- 
tion series of the drop shape under shear. To derive Taylor's Eq.(2.2) it is sufficient to represent 
the drops by their spherical equilibrium shape. Aiming to improve the constitutive equation, 
Schowalter, Chaffey & Brenner (1968) took into account deformations of drops to first order 
in C ■ The refined analysis did not affect the off-diagonal elements of the constitutive equation, 
i.e. Taylor's Eq.(2.2) for the viscosity, but it gave the (unequal) normal stresses to order (9(cy). 
Another limit, where exact results can be obtained, is the limit of large viscosity ratios oo 
(Frankel & Acrivos 1970; Rallison 1980). To clarify the physical significance of the different 
limits we want to give a brief qualitative description of the behavior of a suspended drop under 
shear, based on work by pldroyd (1953 ) and Rallison (198C). 

In a quiescent matrix fluid of viscosity r|t , a single weakly deformed drop relaxes exponen- 
tially into its spherical equilibrium shape; i.e., defining dimensionless deformation by D : = 
{a~ b)/{a + b) with a and b the major and minor axis of the elongated drop, one has for a 
small initial deformation Do, 



D=Doe 



-r/xi 



The characteristic relaxation time ( Pldroyd 1953| ) 

^cR (2A. + 3)(19A.- 



16) 



40(A.+ 1) 



(2.4) 



(2.5) 



also characterizes the macroscopic stress relaxation in an unstrained region of a dilute emulsion. 
At cOTi ~ 1 one observes the characteristic relaxation mode in the frequency dependent moduli. 
The relaxation time diverges for X/o °° since it takes longer for a weak surface tension to 
drive a viscous drop back to equilibrium. What happens if the matrix is steadily sheared at shear 
rate y? For yci ^ 1, the flow induced in the drop by the external driving is weak compared to the 
internal relaxation dynamics and the equilibrium state is only slightly disturbed, i.e., the drop is 
only weakly deformed. Similarly, for large viscosity ratio X, the elongation of the drop becomes 
very slow compared to vorticity, and hence again very small in the steady state, even if Tiy is not 
small. Technically, this is due to the asymptotic proportionality to X^ ' of the shear rate within the 
drop. In both limits of weak deformation, the time Ti also controls the orientation of the major 
axis of the drop with respect to the flow according to 



n 

4 



1 



arctan (Tiy) . 



(2.6) 



Eq.(2.4) and Eq.(2.6) both can be used to determine the surface tension a from observations of 
single drops under a microscope. In passing, we note that the classical method based on the result 
obtained by Taylor ( 1 934 ) for the steady-state deformation, can only be used if X is not too large, 
whereas Eq.(2.4) and Eq.(2.6) are more general. 

The exact calculations mentioned so far became feasible because (and are applicable if) devi- 
ations of the drop from its spherical equilibrium shape are small. On the other hand, if neither 
the capillary number (or Tiy) is small nor the viscosity ratio is large, i.e.. 



C>1 



and 



X< 1 



(2.7) 



drops can be strongly deformed by the symmetric part of the flow field. Experiments with single 
drops by prace (1982 ) and others have shown that this eventually leads to drop break-up if X 
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Figure 1. Schematic representation of the hmits (hatched) where Eqs.(2.13), (2.181 for the viscosity r| 



of an emulsion give rigorous results. In the shaded region they predict r| to be practically independent of 
capillary number. The frame of the box is meant to comprise the whole C-X parameter plane from zero to 
infinity (small Reynolds number understood). The curv ed solid line is a sketch of the break-up curve for 
steadily sheared isolated Newtonian drops according to Grace (1982). Dashed lines indicate schematically 
our viscosity measurements (see Section m. 



is smaller than some critical viscosity ratio. No general rigorous result for the viscosity of an 
emulsion is known in this regime, where arbitrary drop deformations and break-up may occur. 
Below we will also be interested in such cases, where the conditions required for the rigorous 
calculations are not fulfilled. 



2.2. Second-order theory 

For the following discussion we introduce some additional notation. The rate-of-strain tensor 
Cij and the vorticity tensor CO,y are defined as symmetric and antisymmetric parts of the velocity 
gradient djVj. In particular, for a steady simple shear flow v,- ~ jx2?)ii, and 



djVi = cij + = y(5,-i 8^'2 + 5;25/i )/2 + y(5n 8^'2 - 5,-25/i )/2 . 



(2.8) 



Th e components p,/ of the stress tensor in the shear plane {i,j e {1,2}) as obtained for finite X 
by Schowalter et al. (196^ read 



19X+16 

Pij = l^retj -^c<^ll—T—-rrjT——T'Deij 



-ri,.(j)Ti'f- 



25?i2+4a + 4 



\2X + 3)iX+l)^^" ' 14(2X + 3)(X+1)2 

As usual, the material derivative has been defined by 

'Dcij 3, + VkdkCij + (OikCkj + CO jkCii, , 



5/;+0(Yc2). (2.9) 



(2.10) 



where summation over repeated indices is implied, and the first two terms vanish for steady shear 



flow. Note that since ®e,j is diagonal, Eq.(2.9) implies r| = rij- + (9(c ), and hence Eq.(2.2) 



remains valid to first order in C as we mentioned already. 

Can we extrapolate the exact second order result Eq.(2.9) for the stress tensor to arbitrary C 
and X by using the constraints provided by general invariance arguments? For example, since 
the shear stress has to change sign if the direction of the shear strain is inverted whereas the 
normal stresses do not, the shear stress and the normal stresses have to be odd/even functions of 



Y, respectively. From this observation we could have foreseen that Eq.(2.2) cannot be improved 
by calculating the next order in y, i.e., by considering droplet deformation to lowest order More 
important are Gahlean invariance and invariance under transformations to rotating coordinate 
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frames, which give rise to the material derivative introduced above. Applying the operator (1 + 
Xi©) to Eq.(2.9) adds to the right hand side of the equation a term 2r[Tli'Deij plus a term of 
order 0{jC^), so that one obtains (in the shear plane) 

25A.2+41?^ + 4 



Pij + ZiDpij = 2ri7-(e,7 + T2®e;;) +ri,.(|)Ti'y^- 



(2.11) 



14(2?i + 3)(X+l)2 " 

As another short-hand notation we have introduced a second characteristic time T2, which to the 
present level of accuracy in (|) is given by 

19A.+ 16 



12/11 = 1 - (j) 



(2A. + 3)(2A. + 2) 



(2.12) 



It sets t he time scale fo r strain relaxation in an u nstressed region and was named retardation 
time by [Oldroyd (1953 ). Frankel & Acrivos (197C ) realized that up to the partly unknown terms 
of order 0{yc'^) on the right-hand si de, Eq.(2.1 1 ) belongs to a class of possible viscoelastic 
equations of state already discussed by Oldroyd (1958 ). Hence, setting 0{yc^) = on the right- 
hand side of Eq.( 2.11 ), we can d efine a (minimal) model viscoelastic fluid th at beh aves identical 
to the emulsion described by Eq.( 2^ ) for small shear rates. In contrast to Eq.( 2.11 ), the truncated 
formula for the viscosity 



ri=riT 



1 +1:1 



1+(1) 



5A- + 2 

2A. + 2" 



+ 1+^ 



2A. + 3 



(^lY)' 



(2.13) 



thus obtained has a manifestly non-perturbative form. However, no phenomenological parame- 
ters had to be introduced. Note that Eq.(2.13) comprises both exactly known limits: C ^ for 
fixed X, and A, ^ 00 for arbitrary C ■ Obviously, Frankel & Acrivos (1970| ) have forgotten a term 
—25r[c<^ejj /2 in their Eq.(3.6) for p,,- in the limit X ^ °o for fixed C- If the latter is included, 
Eq.(2.13) is also in accord with their (9 (A,^^)— analysis. Moreover, Eq.(2.13) has the proper lim- 
iting behavior for A- = 1, a ^ 0, i.e. C ^ °°, which is an extreme case of Eq.(2.7). Since we 
assume equal densities for the two phases, the two-phase fluid actually reduces to a one-phase 
fluid in this degenerate case, and the viscosity is simply r|c, independent of morphology. For 
illustratio n, the rigorous limits of Eq.( 2.13 ) in the C — A plane are depicted graphically in Fig. |l]. 
Following Grace (1982 ), a qualitative break-up curve for single drops under steady shear is also 
sketched. In summary, Eq.( 2.13 ) is correct for arbitrary A if C ^ 0, and for ar bitrary C if A — > 0°, 
and for small and large C if A = 1. Therefore, one can expect that Eq.(2.13) works reasonably 
well within a large parameter range (small Reynolds number understood). This is further sup- 
ported by the observation that the error made in going from Eq.(2.9) to Eq.(2.13) rather concerns 
the shape of the droplet than its extension (it consists in truncating a perturbation series in shape 
parametrisation). The final result, though sensitive to the latter, is probably less sensitive to the 
former. Nevertheless, one would not be surprised to see deviations from Eq.(2.13) when drops 
become extremely elongated. Finally, due to changes in morphology by break-up and coales- 
cence, the avarage drop size R may change. Observe, however, that for most viscosity ratios (A 
not close to unity), Eq.( 2. 1 3 ) is practically independent of capillary number (and thus of R) when 
break-up might be expected according Grace (1982[ ) and others. As an analytic function that is 
physically known to be bounded from above and from below (the latter at least by the viscosity 
of a stratified two-phase fluid depicted in Fig. H), r|(A, c) has to have vanishing slope in the 
C— direction for large C- According to Eq.( 2.13| ), r|(A, c) is almost independent of C for 

40(A+1) 



C 



(2.14) 



%/3(2A + 3)(19A+16) ' 

where C* is the turning point in the dilute limit, determined by Tiy = 1/ \/3. For finite volume 
fractionsTi from Eq. (2. 16)replaces Xi. Hence, for C C*, Eq.(2.13) and its extension to higher 
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Figure 2. Self-consistent mean-field description of volume-fraction effects in disordered emulsions and 
suspensions. Tlie fluid surrounding a test droplet is assumed to have the viscosity r\c of the continuous 
phase/the viscosity T] of the whole emulsion, within/outside the "free-volume-sphere" of radius R/(^^/^. 



volume fractions derived below in Eq.(2.18) are practically independent of drop deformation 



and morphology. For finite volume fractions a rough interpolation for C* derived from Eg. ( 2. 18 ) 
i s giv en by 2.4/(5(1 +(|)) +4A,). In a large part of parameter space we thus expect Eqs.( 2.13 ), 
( 2.18 ) to be applicable to monotonic shear histories with R given by the average initial radius of 
the droplets. For non-monotonic shear histories, there can of course be hysteresis effects in ri(Y) 
that result from morphological changes for C ^ C*. These can only be avoided by substituting 
for R the radius corresponding to the actual average drop size at the applied shear rate j. 

2.3. Finite volume fractions 

A general limitation of the equations discussed so far, is the restriction to small volume frac- 
tions. Above, we have implicitly assumed that second order effects from drop interactions are 
small compared to second order effects from drop deformation. Any direct (coalescence) and 
indirect (hydrodynamic) interactions of droplets have been neglected in the derivation. Hydro- 
dynamic interactions can approximately be taken into account by various types of cell models. 



Recently Palierne (1990) proposed a self-consistent method analogous to the Clausius-Mossotti 



or Lorentz-sphere method of electrostatics. For the case of a disor dered spatial distribution of 
drops his results reduce to those already obtained by Oldroyd (1953 ). Oldroyd artificially divides 
the volume around a droplet of viscosity r|^/ = Xr[c into an interior "free volume" with a viscosity 
Tj^ of the bare continuous phase and an exterior part with the viscosity r| of the whole emulsion 
(see Fig. ||). According to this scheme, an improved version of Eq.(2^) should be ( Oldroyd 1953 ) 

5 + 3(Tir/Tl.-l) 



5-2(iir/il.-l) 



(2.15) 



This equation predicts a larger viscosity than its truncation to first order in (|), Eq.(2.2). Both are 
shown as dot-dashed lines in Fig. ^. Eq.(2.15) is qualitatively superior to Eq.(2.2). We note, how- 
ever, that the limit A, ^ oo deviates in second order in (|) from the result obtained for suspensions 
by Batchelor & Green (1972). Eq.( |2.15| ) and likewise all of the following equations containing 
quantities fj^, Ti, T2 are only rigorous to first order in (|). 

Th e same reasoni ng as to the viscosity applies to the characteristic times Ti and T2 which now 
read ( |01droyd 19531 ) 

r|,./? [19A. + 16] [2A. + 3 - 2(^(A. - 1 )] 



^2 



a 40(A,+ l)-8(j)(5A. + 2) 

^cR [19A.+ 16][2A. + 3 + 3(j)(A.-l)] 
~o 40(A,+ 1) + 12(|)(5A,-F2) 



(2.16) 



(2.17) 




Figure 3. Comparison of different mixing rules for emulsions with viscosity ratio X = 3 (chosen arbi- 
trarily). The dot-dashed and dashed straight lines pertain to dilute emulsions described by the extrap- 
olation formula Eq.(2.12), which reduces to Tayl or's fo rmula Eq.([2.2|) for small capillary numbers. The 
corresponding curved lines are obtained from Eq.( [2.18[ ) where interactions of the droplets are taken into 
account in a mean-field a pproximation. The curved dotte d lines are the predictions of the cell model by 
Choi & Schowalter (197^ ) for small shear rates, Eq.( [2.19| ). The curved solid line represents the viscosity 
T] = r|c/[(|)+ (1 — (fj^^'j of a two-phase stratified fluid and is a lower bound for any viscosity mixing rule. 



Their ratio 12/^1 is still given by Eg. (2. 12). Finally, Eg. (2. 13) becomes 



ri =ri7- 



1 + (tiY)^ 
ri,. /2X + 2- 



-3(j)(A. + 2/5) 2A. + 3 + 3(j)(A.-l) 



i + \2X + 2-2(i?{X + 2/5) 2X + 3 -2(j)(>i- T 



(^ly)^ 



(2.18) 



which to our knowledge has not been given before, and is one of our main results. (For a graphical 
representation see Fig. Q.) In the limit Tiy ^ it reduces to Eg.(2.15), whereas for Tiy ^ cxj only 
the second term in parentheses contributes and the curved dashed lines in Fig. are obtained. 



From the foregoing discussion one should expect Eg.( 2. 1 8 ) to be applicable within a large range 
of shear rates, viscosities, and volume fractions. 

Finally, we note that more cumbersome expressions for f\T,ti and T2 have been derived within 
a different cell model by [Choi & Schowalter (1975| ). Here we only guote their expression for fjj-. 



Tlr c&s J 



2[(5X)-5(?i-lH^/3] 



4(A. + 1 ) 5 (5A. + 2)(^ + 42A.(^5/3 _ 5 (5A. - 2)(j)V3 + 4(A. - 1 )(j) 'O/^ 



(2.19) 



which is also r epresented graphically by the dotted lines in Fig. ||. Since our data favor Eg. ( 2. 15 ) 
over Eg. ( 2. 19 ), and similar observations have been made by others before (see Section |), we 
will not pursue this alternative approach further in the present contribution. 



2.4. A similarity rule 

It is interesting to observe that if morphology is conserved (drop size R independent of shear rate) 
for C ^ C*, our Eg. (2. 18) for the nonlinear shear viscosity is closely related to the expressions 
for the freguency dependent complex viscosity ri*(co) = ri'(co) +iri"(co) of an emulsion of two 
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incompressible Newtonian liquids as derived by Dldroyd (1953 ) 



. 1+T2ia) 
r|* (0 =117- , , - . 



ri (co) =ri7-- 



l + (Xl(0)2 



For convenience, we also give the corresponding viscoelastic shear modulus G* (co) = icorj* (co) : 
G'(co) + iG"(co), 



- mTTTlT^ ' G"(co) - coiiry— ^ 



l + lTlCO 



(2.21) 



Obviously, the shear-rate dependent viscosity of Eg. ( 2. 18 ) is obtained from the real part of the 
complex frequency dependent viscosity r|* (co) by substituting y for co. 



ri(Y)~ri'(co) 
In the same way, the first normal stress difference 

Pn-P22 = -2Y-(l)r|,.TiY- 



19A.+ 16 



(2A, + 3)(2X + 2) 

fromEq.(§!|) is obtained to leading order in (|) and y from 2cor|"(co), i.e., 

;?ii-/?22^2cori"(co) . 



(2.22) 



(2.23) 



(2.24) 



Reverting the line of reasoning pursued so far, we can conclude that to leading order in C and/or 
"kr^ the weak deformation limit of the constitutive equation of emulsions is obtained from the 
linear viscoelastic spectra G'(co), G"(co). Further, this identification can possibly be extended (at 
least approximately) into regions of the C — A- plane where the critical capillary number for drop 
breakup is somewhat larger than C* of Eq.(2.14). In hindsight, it is not surprising that in the 
case of weakly deformed drops the frequency dependent viscosity and the steady shear viscosity 
are related. Note that under steady shear, drops undergo oscillatory deformations at a frequencey 
2co = Y i f obse rved from a co-rotating frame turning with vorticity at a frequency CO = y/2. If we 
take Eq.( 2.24 ) seriously beyond the rigorously known li mit, w e obta in an interesting prediction 
for the first normal stress difference. In contrast to Eq.( |2.23| ), Eq.( 2.24 ) implies that the first 
normal stress difference saturates at a finite value AQ(Sf<3 /R\Tk + 3 + 2(|)(1 — X)Y for high shear 
rates. Thus, although the initial slope of the first normal stress difference with Y increases with X, 
its limit for large Y decreases with X. 

Finally, we remark that ba sed o n qualitative theoretical arguments, the similarity relation con- 
tained in Eq.( 2.22 ) and Eq.( Z.24 ) has recently been proposed also for polymer melts (Renardy 
1997). Usually, in polymer physics a slightly different rela tion is considered; namely a similarity 
between r|(Y) and |r|*(co)|, also known as Cox-Merz rule (Cox & Merz 1958). In our case, since 
ri"/ri' = G'/G" = 0((|)), we can write 



il(Y)^|Tl*M| + 0(^2 



(2.25) 



Under the conditions mentioned at the beginni ng o f t his se ction, the usual Cox-Merz rule is 
fulfilled to first order in (|) for an emulsion. Eqs.(2.22), (2.24) are interesting from the theoretical 
point of view, because they suggest a similarity of two a priory rather different quantities. The 
results of this section also can be of practical use, since they suggest that two different methods 
may be applied to measure a quantity of interest. 

2.5. Non-Newtonian constituents 
Generalization of the above theoretical discussion to the case of non-Newtonian constituents is 



not straightforward. Indeed, as Dldroyd (1953) already knew, his linear-response results quoted 
in Eq.( 2.2C| ) and Eq.(2.21) are readily generalized to viscoelastic constituents by replacing the 
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Figure 4. Eq.( 2.18 ) normalized to r|c as a function of viscosity ratio X and capillary number C ■ The 
volume fraction of the dispersed phase is chosen to be (|) = 0.3. 



viscosities rj^/ ^ in the expression for r|* (or G*) by complex viscosities r|*y ^.(co) ( ^alierne 1990 ). 
As a consequ ence, the decompositions of r|* and G* in real and imaginary parts are no longer 
those of Eqs.( 2.20 ) and ( 2.21 ), and r|', r|", G', G" are given by more cumbersome expressions. 
For the steady-state viscosity, on the other hand, one has to deal with a non-homogeneous vis- 
cosity even within homogeneous regions of the emulsion, since the strain rate itself is non- 
homogeneous and the viscosities are strain rate dependent. We do not attempt to solve this prob- 
lem here, nor do we try to account for elasticity in the nonlinear case. Yet, it is an intriguing 
question, whether the similarity rule Eq.(2.22) can be generalized to the case of non-Newtonian 
constituents if the constituents themselves obey the Cox-Merz rule (what many polymer melts 
and solutions do). If both constituents have similar phase angles 9 = arctanG"/ G' the generaUzed 
viscosity ratio 



(2.26) 



11.* |r|,*| 

that enters the expressions for rj* and G*, transforms approximately to T|f/(y)/r|t(y) by applying 
the Cox-Merz rule. Therefore, in this particular example, Eq.(2.18) supports the expectation 



that the generalization may work at least approximately. If, on the other hand, the phase angles 
of the constituents behave very differently, the answer is less obvious. This problem has been 
investigated experimentally and is further discussed in Section ^ 

In any case, the generalization can only work if the representation of the emulsion by a simple 
shear-rate dependent viscosity ratio T|f/(y)/r|t (Y), with ythe external shear rate, is justified. In 
the remainder of this section we construct an argument that allows us to estimate the effective 
shear-rate dependent viscosity ratio that should replace A- in Eq.( 2.18 ). We take into account 
the deviation of the strain rate from the externally imposed flow only within the drops, because 
outside the drops the discrepancy is always small. Inside a drop, the strain rate can be small even 
for high external shear rates if the viscosity ratio A, = r|^//r|e is large. Since we are looking for 
an effective viscosity T\d{i) for the whole drop to replace the viscosity r\ci at small shear rate, 
we replace the non-Newtonian drop of non-homogeneous viscosity by an effective pseudo- 
Newtonian drop of homogeneous but shear-rate dependent viscosity. A possible ansatz for fj j is 
obtained by requiring that the total energy dissipated within the drop remains constant upon this 
substitution. Hence, we have 



dV pijgij 



IxXd JdV gijgij , 



(2.27) 
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where gij and gij denote the rate-of-strain fields in the real drop and in the corresponding model 
drop of effective viscosity %t, respectively. They both depend on the position within the drop, 
whereas the average strain rate 



iis=y JdV gijgij 



(2.28) 



that enters the right-hand side of Eg. ( [2. 27 ) does not. Since the strain field gij within a drop 
of homogeneous viscosity is unique for a given system in a given flow field, x\d itself can be 
expressed as a function f\d{%f() of the average strain rate. Here, we approximate this functional 
dependence by the strain rate dependence T\d{i) of the viscosity of the dispersed fluid. Further, 
neglecting drop deformation we calculate Yeff from the velocity field within a spherical drop 



(Bartok & Mason 1958) and obtain 



ii4V2y/(X+1) 



Mi) 



t1c(y) 



(2.29) 



A different prefactor (a/7 in place of V2) in the expression for the effective strain rate Yeff was 
obtained by ie Bruijn (1989) using instead of the average in Eq.(2.28) the maximum norm of 
gij. To obtain the correction to Eq.( [2.18| ) due to Eq.(Z.29) in the case of non-Newtonian con- 
stituents, the implicit equation for A, has to be solved for given functions T\c{i) and 11^/(7). For 
shear thinning constituents, Eq.( [2.29 ) implies a tendency of Eq.(Z.18) to overestimate r|(Y) if y 
and "k are large. In the actual case of interest, for the constituents that were used in the exper- 
iments discussed in Section ^ the viscosity ratio A- and taken at the external shear rate 
Y) varies almost by a factor of 10. However, the corrections discussed in this Section only be- 
come important for high shear rates, where the constituents are shear thinning. In this regime, 
the viscosity ratio (viscosities taken at the external shear rate) only varies between 1/2 and 2, 
and hence the correctio ns exp ected from Eq.( 2.29 ) are at best marginally significant at the level 
of accuracy of both Eq.( 2.29| ) and the present measurements. Therefore, a representation of the 
drops by pseudo-Newtonian drops of homogeneous but shear-rate dependent viscosity is most 
probably not a problem for the measurements presented in the following section. The question 
as to a generalization of the similarity rule Eq.( |2.22| ) to non-Newtonian constituents seems well 
defined. 



3. Experiment 

3.1. Materials and methods 

The experimental investigation deals with a phase separated aqueous solution containing a polysac- 
charide (alginate) and a protein (caseinate). This type of solutions are currently used in the 
food industry. The methods for characterizing the individual polymers in solution are in general 
known, especially when dealing with non-gelling solutions where composition and temperature 
are the only relevant parameters. The polymers are water soluble. When the two biopolymers 
in solution are mixed, a miscibility region appears in the low concentrations range and phase 
separation at higher concentrations. The binodal and the tie lines of the phase diagram can then 
be established by measuring the composition of each phase at a fixed temperature. In general, 
the rheological behavior of phase separated systems is difficult to investigate, and a suitable pro- 
cedure is not fully established. In some cases, two-phase solutions macroscopically separate by 
gravity within a short period of time, but in some other cases (such as ours) they remain stable 
for hours or days without appearance of any visible interface. These "emulsion type" solutions 
have no added surfactant. Following approaches developed for immiscible blends, one may try 
to characterize the partially separated solution as an effective emulsion if the coarsening is slow 
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enough. In order to establish a comparison between phase separated solutions and emulsions, it 
is necessary to know 

• the volume fraction of the phases, 

• their shear-rate dependent viscosities (flow curves), 

• their viscoelastic spectra, 

• the interfacial tension a between the phases, 

• the average radius of R the drops 

Only the ratio R/<3 enters rheological equations. Knowledge of either 7? or a allows the other 
quantity to be inferred from rheological measurements. 

A difficulty when working with phase separated solutions, as opposed to immiscible polymer 
melts, arises from the fact that each phase is itself a mixture (and not a pure liquid) and there- 
fore the rheology of the phase depends on its particular composition. If one wishes to minimize 
the number of parameters, it is important to keep the composition of the phases constant upon 
changing the volume fractions. This can be achieved by working along a tie line of the phase di- 
agram. And this is precisely the procedure that we followed. The polymers were first dissolved, 
then a large quantity of the ternary mixture was prepared (350 ml) and was centrifuged. The two 
phases were then collected separately. Both pure phases were found to be viscoelastic and to 
exhibit shear thinning behavior, which is especially pronounced for the alginate rich phase with 
r|(10^' s^')/r|(10^ s^') w 20, while the caseinate rich phase is almost Newtonian below 10^ 
s^ ' . The viscosity of the alginate rich phase is higher than that of the caseinate rich phase for 
shear rates below 2 • 10^ s^' and lower for higher shear rates. We checked that both phases obey 
the usual Cox-Merz rule in the whole range of applied shear rates. 

By mixing various amounts of each phase, the volume fraction of the dispersed phase was 
varied between 10 % and 30 % while the composition of each phase was kept constant. In partic- 
ular, the temperature was kept constant and equal to the centrifugation temperature in order to 
avoid redissolution of the constituents. To prepare the emulsion, the required quantities of each 
phase were mixed in a vial and gently shaken. Then the mixture was poured on the plate of the 
rheometer (AR 1000 from TA Instruments fitted with a cone and plate geometry 6 cm/2°) and 
a constant shear rate was applied. The apparent viscosity for a particular shear rate was then 
recorded versus time until it reached a stable value. By shearing at a fixed shear rate, one may 
expect to create a steady size distribution of droplets, with a shear rate dependent average size. 
After each shear experiment a complete dynamic spectrum was performed. In this way, shear 
rates rangi ng between 3 ■ 1 0~^ s^' and 10^ s^^ were applied. The analysis of each spectrum ac- 



cording to ^alierne (1990 ) allowed us to derive by curve fitting the average drop radius R at the 
corresponding shear rate. 

More technical details about the experimental investigation along with more experimental 
results will be presented elsewhere. Here, we concentrate on the analysis of those aspects of the 
rheological measurements pertinent to the theoretical discussion in Section |[ 

3.2. Results and discussion 

In this section we present our experimental observations and address the questions posed at 
the end of the introduction. Before we present our own data we want to comment briefly on 



related data recently obtained for polymer melts by arizzuti & Buonocore (1998). These authors 



measured the shear-rate dependent viscosity of binary polymer melts and compared them to the 



low volume fraction limit of Eq.(2.18), i.e. Eq.(2.13), and to (a truncated form of) results of 
Choi & Schowalter (1975). They reported much better agreement with Eq.( [2.13| ) than with the 
truncated series from |Choi & Schowalter (1975 ). Comparison with the full expressions of Choi 
& Schowalter (1975) would have made the disagreement even worse (cf. Fig |). The average 
radius R that enters the equation, was determined independently for each shear rate applied. 
The constituents where moderately non-Newtonian polymer melts, the viscosity ratio varying 
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Figure 5. The nonlinear shear viscosity ri(Y) (opaque squares) and the real part ri'(co) of the dynamic 
viscosity r|*(co) (lines) of a dropl et ph ase of a mixture of weakly viscoelastic polymer solutions (algi- 
nate/caseinate). Also shown is Eq.(2.18) for the viscosity of an emulsion of Newtonian constituents eval- 



uated for the actually non-Newtonian viscosities of the constituting phases with the drop size obtained 
from the spectra (open triangles). Due to the scatter in the dynamic viscosity at low frequencies th ere is an 
uncertainty in the average drop size, resulting in corresponding error bars (multiple points) for Eq.(2.18). 



between A, « 0.3 ... 3 over the range of shear rat es app lied. He nce, th ese experiments, are located 
in the interesting parameter range, where Eqs.( 2.18 ) and Eg. ( 2. 13 ) for r|(Y) are expected to be 
sensitive to drop deformation and break-up. Surprisingly, the results show that they describe the 
data very well over the whole range of shear rates although one would not necessarily expect 
average drop deformation to be very small. Unfortunately, drop sizes have not been reported by 
the authors, so conclusions concerning the location in the C — A, parameter plane and the validity 
of the similarity rule Eq.(2.22) cannot be drawn. Also the question, whether Eq.(2.18) holds for 
small viscosity ratios X <C 0.3, cannot be answered. 

Our own measurements were located in about the same A,— range. As we noted in the preceding 
section, only the ratio /?/o enters rheological equations and kn owledge of either R or a allows the 
other quantity to be inferred from rheological measurements. Ding & Pacek (1999 ) determined 
the interfacial tension of the alginate/caseinate system used in our experiments by observing 
drop relaxation under a microscope and analyzing the data according to Section ^. They found 
a ~ 10^^ N/m. Using this, we obt ained an averag e drop size R ~ 10^^ m from the measured 
spectra G'((b), G"{(£i) according to Pali erne (199C ) for the experiments reported in Fig. ^j. By 
the method based on Palierne (199C), we could not detect a decrease in drop size with shear 
rate as expected from the phenomenological phase diagram for single Newtonian drops under 
shear as established by Grace (1982) and others. Thus, the limit of high capillary numbers and 
moderate viscosity ratios (the region above the break-up curve) in Fig. |l] has been accessed 
experimentally. Corresponding locations have been indicated qualitatively in the figure by dashed 
lines. With the drop size being constant, one can try to test the proposed similarity rule Eq.(2.22). 
By identifying the axis for frequency co and shear rate y, data for the real part r|'(a)) of the 
frequency dependent dynamic viscosity r|* (co) are compared to data for the shear-rate dependent 
viscosity r|(Y) in Fig. |[ The emulsion containing 30% of the alginate rich phase and 70% of 
the caseinate rich phase has been prepared at room temperature as described in the preceding 
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section. Steady shear rates ranging between 3 • 10^^ s^' and 10^ s^' corresponding to capillary 
numbers C ~ 6 • 10^^ ... 10^ have been applied. The shear viscosity (opaque squares) is reported 
in the figure for each of these individual measurements. The multiple data sets for rj' (lines) 
taken each between two successive steady shear measurements, superimpose fairly well; i.e. 
the spectra appear to be remarkably independent of the preceding steady shear rate. The good 
coincidence of r|'(a)) and r|(y) in Fig. ^ show that the data obey the proposed similarity rule 
Eq.(2.22) over a large range of shear rates. The agreement near 2 • 10^ s^' is a consequence 
of the proximity to the trivial limit X—\,C—°°- Nevertheless, the data provide strong evidence 
that Eq.( 2.22 ) is an excellent approximation for a large range of viscosity ratios and capillary 
numbers. Simi lar results (not shown) have been obtained for other volume fractions. Comparison 
with Eq.( 2.18 ) represented by the open triangles in Fig. ||, on the other hand, is less successful at 
large shear rates, although it is still not too far off for a theoretical curve without any adjustable 
parameter. A discrepancy had to be expected as a consequence of the non-Newtonian char acter 
of the constitu ents a t large shear rates, which is definitely not taken into account in Eq.( 2.18 ). For 
the plot of Eq.( 2.18 ) in Fig. ^ we merely substituted ^diy) /^cH) taken at the external shear rate 
Y for X. The average drop size R was obtained from fitting the viscoelastic moduli. The scatter 
in the dynamic viscosity data gives rise to an uncertainty in R, which is reflected by the multiple 
open triangles at low shear rates. It seems that the similarity rule Eq.(2.22) is more general than 
Eq.(2.18), i.e., it still holds for rather viscoelastic constituents (that obey the usual Cox-Merz 
rule), where the latter fails. This relation certainly deserves further investigation with different 
materials and methods. 

In summary, we have succeeded in establishing an analogy first between a partially phase- 
separated polymer solution and an emulsion, and further between the viscoelastic spectrum of 
the system and its nonlinear shear viscosity even in the case of (moderately) non-Newtonian 
constituents. 
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